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Abstract. I review the topic of relic neutrino asymmetry generation through active- 
sterile neutrino and antineutrino oscillations in the early universe. Applications to (i) 
CO . the suppression of sterile neutrino production, and (ii) the primordial Helium abun- 

dance are briefly presented. 
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Reasonably large relic neutrino asymmetries will be generated by active-sterile 
neutrino and antineutrino oscillations in the early universe provided that certain 
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mild conditions are met [1]. "Reasonably large" in this context means that the 
neutrino asymmetry, defined by 

' T _ n Ua - n 

where a = e, /i, r and the n's are number densities, can have a final value as large 
as about 3/8 [2,3]. The mild conditions are that the vacuum mixing angle 6q should 
be small and that the Am 2 should be negative. The flavour and mass eigenstates 
are related by 

\v a ) = cos9 \v a ) + sm9 \v b ) (2) 
\u 8 ) = -sin^oK) + cos9 \u b ), 

so Am 2 < means that the predominantly sterile mass eigenstate is lighter than the 
predominantly active mass eigenstate. We focus on that epoch of the early universe 
subsequent to the disappearance of a significant muon/antimuon component to the 
plasma, and ending with the period of Big Bang Nucleosynthesis (BBN). 

One reason the generation of neutrino or lepton asymmetries is interesting is 
evident from the effective matter potential, 
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(3) 



where p is the neutrino momentum, and 
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The function a is just the generalisation of the usual Wolfenstein effective potential 
[4] pertinent to the considered epoch of the early universe. Gf is the Fermi constant, 
T is the temperature and the "effective asymmetry" is given by 

L^ = L Va +L Ve + L Vll + L Vr +ri, (6) 

where rj ~ 1CT 10 — 10~ 9 is related to the baryon and electron asymmetries required to 
produce the universe we observe. When lepton asymmetries are large, the effective 
potential is also large, leading to small matter-affected mixing angles and thus 
also to the suppression of sterile neutrino production. The asymmetry generation 
effect must be taken into account when exploring the implications for BBN of 
neutrino scenarios involving light sterile degrees of freedom. In particular, <-> 
v s oscillation parameters motivated by the atmospheric neutrino deficit are far 
from being necessarily ruled out by BBN [1,2,5,6]. The function b is the leading 
finite-T correction to the Wolfenstein term [7]. The constant mw is the H^-boson 
mass, while A a is a numerical factor given by A e ~ 17 and A^ T ~ 4.9. At high 
temperatures, b is large enough to suppress sterile neutrino production by itself. 
However, it decreases with temperature as T 6 , so a large asymmetry is eventually 
required to keep sterile neutrino production suppressed at lower temperatures. 
Another reason lepton asymmetries are interesting is conveyed by the reactions, 

v e n <-> e~p, V e p <-> e + n, (7) 

which play an important role in setting the n/p ratio at weak freeze-out just prior 
to BBN. A sufficiently large electron neutrino asymmetry (a few percent) will alter 
this ratio and hence change the predicted Helium abundance [3,8,9]. In particular, 
successful BBN can be achieved with a higher than usual baryon density, provided 
that L Ue is positive and of the correct magnitude. The threatened overproduction 
of Helium is compensated by the preferential conversion of neutrons into protons 
due to the positive L Ve . 

Light sterile neutrinos are themselves of considerable interest independent of their 
possible role in asymmetry generation. For instance, two-fold maximal mixing is 
easily achieved through a pseudo-Dirac structure [10] or through ordinary-mirror 
neutrino mixing [11]. (Mirror neutrinos are strictly speaking not sterile because 
they feel mirror weak interactions. However, they are effectively sterile from the 



point of view of ordinary weak interactions.) Two-fold maximal mixing is of great 
interest in light of the solar and atmospheric neutrino deficits. In addition, it 
is well known that least one light sterile flavour is required if oscillations are to 
simultaneously explain the solar, atmospheric and LSND neutrino data. 

As a final comment about motivation, let me also add that collision and matter 
affected neutrino oscillation dynamics is a fascinating subject in its own right. As 
we will see, neutrino asymmetry generation is driven by a subtle interplay between 
quantal coherence and decoherence. 

Consider for simplicity a two-flavour active-sterile system v a + v s together with 
its antiparticle counterpart. The appropriate dynamical variables are the 1-body 
reduced density matrices for the neutrinos and antineutrinos, 

p = ^(P + P-a), p=±(P + f.<r), (8) 

respectively. The decomposition with respect to the 2x2 identity matrix and the 
Pauli matrices proves to be convenient. All of these quantities are functions of 
neutrino momentum p and time t (or equivalently temperature T). 

The diagonal entries p are appropriately normalised distribution functions for v a 
and v s : 

Paa = l(Po + Pz) 
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Pss = -(P -Pz) = N eq( y ( 9 ) 

where rif — Nfdp and N eg (£) is the Fermi-Dirac (FD) distribution with dimen- 
sionless chemical potential £ = p Ua /T, 

- 2 ^-0 + r <10) 

We have chosen to normalise the neutrino distribution to the FD distribution with 
zero chemical potential. 
The off-diagonal entries, 

Pa. = pl> = \(Px-iPv), (11) 

are coherences. They quantify the amount of quantal or phase coherence enjoyed by 
the neutrinos. These quantities are needed because non-forward neutrino scattering 
off the background plasma decreases quantal coherence. The entries of p have a 
similar interpretation. 

The evolution of p is given by the Quantum Kinetic (or Boltzmann) Equations 
(QKEs) [12], 
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where D = T/2 and Pp = P^rc + P y y with x, y, z being unit vectors in the stated 
directions. The second equation above has an approximate equality sign because the 
righthand side assumes that all background fermions have thermal FD distributions, 
and that v a is distributed in an approximately thermal manner. 

The V x P term is just a re-expression of coherent matter- affected oscillatory 
neutrino evolution. It is equivalent to the usual Schrodinger Equation governing 
matter-affected neutrino oscillations. Note, however, that in the early universe this 
evolution is non-linear because of neutrino scattering off the background neutrinos 
of the same flavour. The vector V is given by 



V = (3x + \z, (13) 



where 
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/5=^— sin 2^ , \ = -—- cos 26 + V eff . (14) 

The non-linear effects enter through V e ff. 

The —DPt term drives collisional quantal decoherence. The decoherence rate is 
equal to half of the collision rate T for a neutrino of momentum p, where 

T^y a G 2 F T 5 -^. (15) 
(P) 

The quantity (p) ~ 3.15T is the average momentum for a FD distribution with zero 
chemical potential, while y e ~ 4 and y^ T ~ 2.9. This expression for T is correct 
provided that thermal equilibrium holds and asymmetries are not very large. The 
—DPt term tries to exponentially damp the coherences P XjV to zero. Since the 
collision rate goes as T 5 , quantal coherence is expunged at high temperatures. In 
typical applications, the V x P term begins to dominate over the —DPt term as 
T approaches a few MeV. 

The dPo/dt equation describes the repopulation of the u a distribution from the 
background plasma heat bath. It is proportional to the total weak collision rate T, 
and the term in square brackets on the righthand side acts to drive the actual v a 
distribution function (P + P z )N eq (0)/2 to FD form N eq (£). 

The neutrino asymmetry L Ua enters these equations through the function a in 
V e jf, and through the chemical potential £ in 8P /dt. For small L Ua , 



In thermal equilibrium the existence of a neutrino asymmetry is equivalent to the 
existence of a nonzero £ (where the chemical potential for antineutrinos is equal 
and opposite that of the neutrinos above the chemical decoupling temperature). 

The antineutrino QKEs are of the same form as the above equations with the 
substitutions L Ua — > — L Va and — > —L^ a \ 

We are primarily interested in the evolution of the asymmetry L Va . This is 
indirectly given by the neutrino and antineutrino QKEs. It is easy enough, how- 
ever, to derive a redundant but useful direct evolution equation for the asymmetry. 
Employing the QKEs and a + s lepton number conservation one obtains 



This equation is numerically useful because, during most of its evolution, L Ua is the 
difference of two large numbers (the neutrino and antineutrino number densities per 
photon). The use of Eq.(17) circumvents this numerical difficulty This equation 
is also a very useful starting point for developing approximate evolution equations 



We now discuss the behaviour of the QKEs. We will suppose that the initial 
asymmetries are very small or zero, and that a negligible fraction of the primordial 
plasma is in the form of sterile states. Also, we consider |Am 2 | values higher 
than about 10 -4 eV 2 , because for lower values asymmetry growth begins when 
decoherence is negligible, and it tends to be oscillatory [14]. 

For sufficiently high temperatures T, neutrino and antineutrino oscillations are 
severely damped or frozen. There are two reasons for this. First, the decoherence 
function D ~ T 5 is large and drives P XjV and P x y to zero, and hence also the 
righthand side of Eq.(17) (Quantum Zeno Effect). Second, the 6-term in V e ff also 
rises as T 5 , so the effective matter mixing angle is very small anyway. 

As T decreases with the expansion of the universe, the oscillations begin to 
unfreeze. Their initial non-trivial evolution is, however, still dominated by non- 
forward scattering. During this phase, a very interesting phenomenon will occur 
provided that the mild conditions discussed at the beginning of this article are met: 
small 6*o and Am 2 < 0. There is a critical temperature T c , roughly given by 



Above T c , the evolution equations drive L Va such as to impel the effective asym- 
metry towards zero (in other words, L Ua evolves from zero such as to cancel 
the rj term in the effective asymmetry). We can call this the asymmetry destruc- 
tion phase, and = is a stable fixed point of the evolution equations. At T c , 
this stable fixed point becomes unstable and runaway positive feedback leads to an 
exponential increase in L Ua during the explosive growth phase. The approximately 
exponential growth is cut off by the non-linear terms in the QKEs after L Va reaches 




(17) 



[2,13]. 



T c ~ (15 -> 18 MeV) cos2fl ( 
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FIGURE 1. Examples of lepton asymmetry growth curves driven by v r <-> v s and the corre- 
sponding antineutrino oscillations. The mixing angle is selected to be sin 2 20 = 10~ 8 . The three 
curves correspond to Am 2 = —0.5, —50, —5000 cV 2 , reading left to right. This figure is taken 
from Rcf. [5]. 

a value which is several orders of magnitude higher than the baryon and electron 
asymmetries. A power law growth phase, as per L Ua ~ T -4 , then sets in. Typically, 
collision dominated evolution gives way to oscillation dominated evolution during 
the power law phase. This is important, because physically there is a "hand over" 
from collisional asymmetry growth to non-linear MSW driven growth during this 
phase. Finally, asymmetry growth stops when L Ua reaches about 0.2 — 0.3. 

Some examples of asymmetry evolution are given in Fig.l, which is taken from 
Ref. [5]. These curves were produced by numerically solving the QKEs together 
with Eq.(17). An analytical understanding of why these curves take the displayed 
form has been given in the literature [1-3,13]. The key idea is to take the adiabatic 
limit of the QKEs. At high T, it is then possible to define precisely what is meant 
by collision dominated evolution, and to derive an approximate dL Ua /dt equation 
which makes the critical behaviour at T c manifest. At low T, when collisions can 
be neglected, the adiabatic limit of the QKEs is exactly the same as the adiabatic 
limit for non-linear MSW evolution. See Refs. [2,3,13] for a complete discussion. 

I will now review two applications, beginning with the suppression of sterile 



neutrino production. Suppose that we wish to solve the atmospheric neutrino 
problem by maximal — > v s oscillation with |Am 2 s | in the range 1CT 3 eV 2 to 1CT 2 
eV 2 . 3 If the v il J rv s subsystem does not mix with any other neutrino, then the sterile 
neutrinos (and antineutrinos) will certainly be brought into thermal equilibrium 
with the rest of the plasma prior to BBN (unless there are large pre-existing lepton 
asymmetries [16]). The resulting increase in the expansion rate of the universe will 
lead to a higher than standard Helium-4 abundance. I will assume for the sake 
of the example that primordial abundance observations cannot tolerate such an 
increase in the expansion rate. (A complete account of the somewhat complicated 
status of primordial abundance observations vis-a-vis BBN is beyond the scope of 
this talk.) 

Small mixing between the v s and, say, a more massive v T can completely change 
the conclusion that v s is brought into thermal equilibrium. The v T jv s oscillation 
parameters, as chosen in the previous sentence, satisfy the requirements for large 
L Vt generation. If the large L Vt does not have its effects cancelled (see below), 
and if it is generated early enough, then the large matter potential consequently 
generated for the <-> v s mode will suppress these oscillations and hence also v s 
production. A numerical calculation is required to properly analyse the outcome, 
because the maximally mixed <-> v s mode always tends to induce — > 0. 
In other words, the large L Vr that is being created by the u T <-> u s mode could 
be compensated by L v creation driven by <-> u s , with the overall effect being 
that the /x-like effective asymmetry is driven to zero, and consequently that v s 
production through the channel is unsuppressed after all. 

Figure 2 shows the outcome of such a numerical calculation [2,5,6]. It is a 
plot in v T jv s oscillation parameter space. The solid lines correspond to |Am 2 J = 
10~ 3 , 10~ 2 ' 5 , 10~ 2 eV 2 in ascending order up the page. Choose your favourite at- 
mospheric Am 2 s . In the parameter region above the relevant solid line, the L Ut 
asymmetry is not cancelled such that — > 0. Below the solid line, it is can- 
celled. The transition is very sharp. So, above the solid line, v s production via 
the — > v s mode is very suppressed. The dot-dashed line refers to v s production 
via the other available mode, v T — > v s . This small angle mode receives very little 
"self-suppression" from its own L Vt creation activity, so its oscillation parameters 
must satisfy a different sort of bound. This bound is actually very similar to the old 
constraints calculated in Ref. [17] with asymmetry generation artifically switched 
off. For illustrative purposes, the dot-dashed line assumes that no more than 0.6 
extra effective neutrino flavours are allowed. 

The second application concerns the direct effect of a v e asymmetry on the neu- 
tron to proton ratio at weak freeze out. The implications of neutrino asymmetry 
generation in multi-flavour scenarios must be studied on a case-by-case basis. Here 
I will consider the scenario of Ref. [18]: there is one sterile flavour, which is used to 
solve the solar neutrino problem via an MSW style solution. The more massive 

3 ) SupcrKamiokande claim that present data disfavour the — ► v s possibility relative to the 
Up — > v T mode [15]. However, the dust has yet to settle on this issue. 




FIGURE 2. Allowed region according to BBN for v T <-> v s oscillation parameters in order to 
suppress maximal <-» v s oscillations via a large L Vt . See text for a complete discussion. This 
figure is taken from Ref. [5]. 



and v T states are maximally mixed, with the mass gap between (z/ e , v s ) and (u^, v T ) 
set in the LSND range. In this scenario, both the — ■> v s and v T — > z/ s modes 
satisfy the conditions for asymmetry generation. Once L„ M and Lj, T have been gen- 
erated, these asymmetries can be reprocessed into L Ve via <-> z/ e oscillations 
[8]. The outcome is insensitive to the v e /v^ T mixing angles for a large range 
of these parameters because the MSW transitions are adiabatic. It is, however, 
somewhat sensitive to Am 2 between (z/ e , u s ) and (z/ M , u T ). Figure 3 shows a plot of 
the change in the effective number of neutrino flavours during BBN as a function 
of this Am 2 parameter [8]. The "effective number of neutrino flavours" is just a 
convenient measure of the change in the Helium mass fraction Y p relative to stan- 
dard BBN, roughly obeying the relation 5Y P ~ 0.0125 . It is affected by both 
expansion rate alterations and L Ue . Figure 3 incorporates both effects. 

There are a few other interesting neutrino models that have been analysed in 
the literature. A particularly interesting case is the Exact Parity or Mirror Matter 
Model [11]. The full analysis is very complicated; please see Ref. [19] for all of the 
details. 

To conclude: If light sterile or mirror neutrinos exist, then they should play a ma- 
jor role in cosmology. The oscillation generated neutrino asymmetry phenomenon 
would be a central feature. Sufficiently large neutrino asymmetries would suppress 
sterile or mirror neutrino production, and an electron-neutrino asymmetry of the 
right magnitude and sign would affect primordial Helium abundance. In particular, 




FIGURE 3. The change in the effective number of neutrino flavours during BBN for a certain 
neutrino scenario, plotted as a function of a Am 2 parameter that was taken to be in the LSND 
range. See the text for a more complete discussion. This figure is taken from Ref. [8] 



a positive L Ve would allow a larger baryon density to exist without the overpro- 
duction of Helium. We note with interest that the recent Boomerang/Maxima 
cosmic microwave background anistropy measurements favour a larger than stan- 
dard baryon density [20]. While the combined solar, atmospheric and LSND data 
require at least one light sterile neutrino if oscillations are to simultaneously re- 
solve all of the anomalies, we especially look forward to future experiments that 
could provide further evidence for light sterile neutrinos: SNO, MiniBooNe, the 
longbaseline experiments, as well as further SuperKamiokande data. 
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